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ABSTRACT 

We report on a new Be/X-ray pulsar binary located in the Wing of the Small Magellanic Cloud 
(SMC). The strong pulsed X-ray source was discovered with the Chandra and XMM-Newton 
X-ray observatories. The X-ray pulse period of 1062 s is consistently determined from both 
Chandra and XMM-Newton observations, revealing one of the slowest rotating X-ray pul- 
sars known in the SMC. The optical counterpart of the X-ray source is the emission-line star 
2dFS 3831. Its B0-0.5(III)e-H spectral type is determined from VLT-FLAMES and 2dF optical 
spectroscopy, establishing the system as a Be/X-ray binary (Be-XRB). The hard X-ray spec- 
trum is well fitted by a power-law with additional thermal and blackbody components, the 
latter reminiscent of persistent Be-XRBs. This system is the first evidence of a recent super- 
nova in the low density surroundings of NGC 602. We detect a shell nebula around 2dFS 3831 
in Ha and [O m] images and conclude that it is most likely a supernova remnant. If it is 
linked to the supernova explosion that created this new X-ray pulsar, its kinematic age of 
(2 - 4) X 10^ yr provides a constraint on the age of the pulsar 
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1 INTRODUCTION 

The Small Magellanic Cloud (SMC) is host to a large population of 
~50 high-mass X-ray binaries (HMXBs), comparable to the num- 
ber known in the Galaxy (e.g. iHaberl & Pietschll2004l : [ Coe et al.l 
l2005l : lReiall20Tlh . However, unlike the Galactic population, all of 
these systems but one are Be/X-ray binary (Be-XRB) systems, with 
the e xception (SMC-Xl) containing a supergiant companion in- 
stead ( IWebsteretai]|l972l) . 

The Wing of the SMC is the region linking the eastern side 
of the SMC to the Bridge. Despite its lower content of gas, dust, 
and stars, it had a major star forming event ~ 1 1 Myr ago and ap - 
pears, as expected, deficient in Be-XRBs jAntoniou et al.ll2010t) . 
Note however that the coverage of the SMC Wing by X-ray obser- 
vations has been sparser than that of the Bar. Each new discovery 
of a Be-XRB in the SMC W ing is therefore particularly noteworthy 
(e.g. lMcGowan et alj|2007h . 

We have conducted Chandra and XMM-Newton observations 
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of the far eastern region of the Wing, concentrating on the star 
forming region NGC 602. The full description of these observa- 
tions will be presented in a forthcoming paper (Oskinova et al., 
in preparation). In this letter, we concentrate on one bright X- 
ray source located near NGC 602 in projection and coinciding 
with the emission-line star 2dFS3831 (R.A. =l''27'" 46.03', Dec. 
= -73°32'56.42", J2000.0), which is revealed as a Be-XRB. Based 
on the results of section 12.21 throughout this paper we will refer 
to the X-ray sou rce as SXP 1062 following the nomenclature of 
I Coe et al.l ( l2005h for SMC X-ray pulsars. 

The stellar population of NGC 602 (which comprises three 
star clusters: NGC602A, B, and C) has been studied in detail. In 
addition to very old stars of 6-8 Gyr, likely the SMC field pop- 
ulation, the NGC 602 clusters contains yo u ng stars of ~4-5 My r 
jWesterlundl 1 1964 iHutchings et"ai] Il99ll : ICignoni et all \200% . 
The significant num ber of young stella r objects detected in the 
NGC 602A cluster jCarlson et all l2007l : iGouliermis etai] l2007h 
suggests that star formation is still ongoing. As for what initially 
triggered the formation of these clusters, the low density environ- 
ment in the Wing of the SMC hints at an additional mechanism 
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acting along with gravitational colla pse. NGC 602 is lo cated at the 
intersection of three H i shells, and iNigra et alj l llOOSh suggested 
that its formation is the result of the in t eractio n of two expand- 
ing shells ~7 Myr ago. ISchmalzl et al. proposed instead 
that star formation was possibly induced by encounters with the 
Large Magellanic Cloud or the Milky Way. NGC 602 is associated 
with SGS-SMCl, the only Ho- supergiant shell known in the SMC. 
The discovery of a Be-XRB in this area (only ~7' to the west of 
NGC 602A) provides important clues to the star formation history 
in the Wing and the evolution of large-scale structures in the inter- 
stellar medium. 

This letter is organized as follows. We describe the X-ray ob- 
servations and analysis in section |2] In section [3] we present the 
spectroscopy of the optical counterpart. We discuss the properties 
and environment of SXP 1062 in section|4]and present our conclu- 
sions in section|5] 



2 X-RAY OBSERVATIONS 

The X-ray data were obtained with the ACIS-I camera on the 
Chandra X-ray Observatory and with the EPIC cameras onboard 
XMM-Newton. The Chandra observations consisted of 1 1 separate 
exposures acquired between 2010-03-31 and 2010-04-29 (effec- 
tive exposure time of 290.7 ks), while 4 separate exposures were 
obtained with XMM-Newton between 2010-03-25 and 2010-04- 
12 (EPIC-pn effective exposure time of 162.5 ks). The data were 
reduced using the most up-to-date versions of the correspond- 
ing data reduction software. The X-ray source CXO J012745.97- 
733256.5 (=SXP 1062) coinciding with the optical emission-line 
star 2dFS3831 is detected in each of these 15 exposures. In the 
Chandra observations, the positional uncertainty (1 cr) of the X- 
ray source is 0.9". The source was also detected during an XMM- 
Newton slew on 2009-11-16. Interestingly, the source was not de- 
tected in the ROSATAll Sky Survey, although at its present luminos- 
ity it would have been bright enough at an expected ROSAT PSPC 
count rate of 0.01 counts s"' . 



2.1 Spectral analysis 

Figure[T] shows the combined background-subtracted Chandra and 
XMM-Newton spectra of SXP 1062. The Chandra spectra were ex- 
tracted from a 12'.' 8 circular region and the background was ex- 
tracted from a concentric annular region of radii 12'.' 8 and 19'.' 2 
which was free of sources. Likewise, the XMM-Newton spectra 
were extracted from a 32" circular region and the background from 
4 adjacent circular regions free of sources with radius 45". 

As a first step the background-subtracted Chandra and XMM- 
Newton spectra of SXP 1062 were simultaneously fitted using a 
simple absorbed power-law model. The best-fit parameters, listed 
in Table[T] imply an absorbed flux in the energy range 0.2-12.0 keV 
of /x= 1.8x10"'^ erg cm"^ s"'. This corresponds to an intrinsic X- 
ray luminosity in this same energy range of Lx=6.9xl0''^ erg s"' 
assu ming a distance mo dulus of 18.7 appropriate for the Wing 
(e.g. ICignoni et al.ll2009l) , or Lx=8.2xl 0^^ erg s~' for the 'stan- 
dard' SMC distance modulus of 18.9 (e.g. lHarries et al.l2003l) . The 
photon index T of ~0 .75 is typical of X-ray pulsar binaries (e.g. 
lYokogawa et al.ll2003l) and a signature of accretion onto a strongly 
magnetized neutron star. 

The simple power-law model suggests some emission excess 
below I keV and at high energies. As a next step, we added a black- 
body component, which statistically improved the fit (see Table[Tll. 
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Figure 1. Chandra ACIS-I and XMM-Newton EPIC-pn spectra of 
SXP 1062. Overplotted is the best-fit model including power-law, thermal, 
and blackbody components. The lower panels show the residuals of the fit. 



The observed X-ray flux is similar to that of the previous model, but 
the X-ray luminosity decreases by ~ 10%. We also tried a model 
for which a thermal plasma component with SMC abundances is 
included instead of a blackbody. Similarly to the blackbody, the 
addition of a thermal component to the power-law model improves 
the fit. In this case, the intrinsic X-ray luminosity is only marginally 
greater by ~ 1% compared to the simple power-law model. 

Finally, the best fit is achieved by adding simultaneously a 
thermal and a blackbody component to the power-law (see Ta- 
ble[T]and Figure[Tll. The flux for this model (PBT in Table [T]! is 
/x= 1.66^025X10"'^ erg cm"^ s"', corresponding to an intrinsic lu- 
minosity Lx =6.3^55gXlO''' erg s"' for a distance modulus of 18.7. 

The column densities are well constrained by spectral fitting, 
albeit the values slightly differ depending on the model. Using the 
best-fit values, we note that the H column densities in Table[T]imply 
Ay = 0.63 - 0.87 mag for the A'h/Av ratio of 1.87 x 10^' cm^^ i.e. 
Eb-v = 0.20 - 0.28 assuming Ry = 3.1. This range is in good 
agreement with the Eb-v = 0.19 value c omputed by c omparing 
the value (B - V) = -0.04 of 2dFS383I ( lMasse\]|2002l) with the 
intrinsic value (S - V)o = -0.23 of a B0-0.5III star l lWegneill 1994 
given the uncertainties on the photometry and spectral type (see 
section[3]for the optical spectroscopy of 2dFS 3831). 



2.2 Timing Analysis 

The photon arrival times were corrected for the solar system 
barycenter. We searched for pulsations in the X-ray light curves 
in the soft and hard energy bands (0.4 keV — 1.5keV, 2.0 keV — 
7.2 keV) and a total energy band (0.4 keV — 7.2keV), using fast 
Fourier transform and light-curve folding techniques as imple- 
mented in the timing analysis software xronos. 

Figure|2]shows the inferred Chandra and XMM-Newton power 
density spectra in the soft and hard energy bands with a clear peak 
at a frequency of 9.4 x 10"'' Hz (P= 1062 s). This coherent X-ray 
pulse period establishes the source as a binary X-ray pulsar. The 
pulse profiles folded with this period are shown in Figure[3] 

The light-curves used to compute the XMM-Newton power 
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Table 1. Best-fit parameters of the spectral models. The code for the different components of each model is: "P" for power-law, "B" for blackbody, and "T" 
for optically thin thermal emission with SMC abundances. 
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Figure 2. XMM-Newton and Chandra power density spectra in the soft and 
hard energy bands for different bin times. 



spectrum in the hard energy band, the XMM-Newton power spec- 
trum in the soft energy band, and the Chandra power spectrum in 
the hard energy band were binned by respectively 5 s, 200 s, and 
200 s. The 5 s bin time for one of these was chosen to make sure that 
a shorter pulse was not missed. No pulse was found in the soft band 
Chandra observations, possibly because the count rate is too low. 
Because of the softer response of the EPIC-pn camera, the pulse in 
the soft band is the most obvious in the light-curves obtained with 
this instrument. 

Note that in addition to the pulse period. X-ray variability was 
also detected at a level of ~ 20% peak-to-peak on a timescale 
of several days in both XMM-Newton and Chandra light-curves. 
Given that this variability does not appear regular and that the time 
sampling is very sparse, we do not analyse these variations further. 



3 OPTICAL SPECTROSCOPY 

From its X-ray position with subarcsecond precision from Chan- 
dra images , the optical coun terpart to SXP 1062 is identified as 
2dFS 3831 tevans et al j2004h. We observed this star with the VLT 
FLAMES instrument jpasguini et all |2002|) on 2010 October 25 
as part of a spectroscopic survey of massive stars in NGC 602 
complementing our X-ray observations. Spectra were obtained in 
the Medusa-fibre mode of FLAMES using the LR0 2 setting of 
the G iraffe spectrograph (3960-4564 A, R=7000, e.g. lEvans et"al] 
l201lh . Five pairs of 1800 s exposures were obtained. The standard 
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Figure 3. Pulse profiles folded with a period of 1062 s in different energy 
bands. 



data processing (bias subtraction, fibre location, summed extrac- 
tions, division by flat-field, wavelength calibration) was done us- 
ing the ESQ Common Pipeline Library FLAMES reduction rou- 
tines (v.2.8.7). Additionally, heliocentric correction and subtrac- 
tion of a median s ky spectrum were performed (for details see e.g. 
lEvans et"al]|201 ih . The spectra from all the exposures were then 
normalized and merged. The 2dF spectra of this star were also re- 
trieved. These two spectra cover the regions from ~3800 to 4900 A 
and -6000 to 7000 A, with re solving powers of 1500 and 2500 re- 
spectively tevans et al.ll2004l) . Figure H shows the VLT FLAMES 
and 2dF spectra of 2dFS 3831. The VLT FLAMES spectrum has 
been smoothed and rebinned to an equivalent resolving power of 
R=4000. 

A weak He n /14542 absorption line is visible in the VLT 
FLAMES spectrum and there is a hint of a weak He ii /14686 
line in the 2dF spectrum, but He ii /i4200 is absent, suggesting a 
spectral type aroun d BO-0.5 following the classification adopted by 
tevans et al.l j2004h . 

Several characteristics attributable to a circumstellar disc in- 
dicate that 2dFS3831 is a classical Be star. Fe ii /i4 179 and Fe ii 
/i423 3 emission lines typical of some early Be stars (e.g. ISlettebakI 
Il982h are visible in the VLT FLAMES spectrum. The Ha emission 
seen in the 2dF spectrum (see the right panel in Fig.|4l( is relatively 
strong with an equivalent width (EW) of -23 A. There is significant 
H/3 emission, plus emission in the core of the other Balmer lines 
(double-peaked in the VLT FLAMES spectrum) and apparent in- 
filling of the He i absorption l ines. Two Micron All Sky Survey 
(2MASS ; ISkrutskie et .^1120061) JHKs photometry also indicates a 
clear infrared excess when compared with expected colors for a 
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Figure 4. Optical VLT FLAMES and 2dF spectra of 2dFS 3831 with key spectral features identified. The VLT FLAMES spectrum is ofl'set vertically. 



BO-0.5 spectral type ( IWegnei|[l994l) . The double-peaked emission 
in the Balmer lines is not due to oversubtraction of nebular features 
present in the median sky spectrum. Also, given the relative weak- 
ness of the nebular emission lines in the mean sky spectrum, the 
emission in the core of the Balmer lines is most likely dominated 
by circumstellar material and not by nebular conta mination. 

The observed magnitude (V= 14.36, B= 14.32: lMassevll2002l) 
and a range of reasonable extinction estimates (including both the 
interstellar and circumstellar compone nts) lead to an absolute mag- 
nitud e consistent with a BO-0.5 giant jWalbomll 19721 : IVacca etal] 
Il996h . We thus determine the spectral type of 2dFS3831 as 
B0-0.5(III)e-l-, where '+' signifies the presence of Fe ii. This 
corresponds to a typical effective temperature T^ff ~ 26 000 K 
and an evolutionar y mass M ~ 15 M© at SMC metallicity (cf. 
iTrundle et alj2007l) . 

We estimated the radial velocity of 2dFS 3831 from the VLT 
FLAMES spectrum by fitting Gaussian profiles to the wings of Hy, 
US, He, He i /14143, and He i /14388. The average of all measure- 
ments is 167 km s^'with a standard deviation of 13 km s"'. This 
is consistent with the me an velocity of massive stars in the SMC 
tevans & HowarthI 20081). a nd also with the gas velocities mea- 
sured by iNigra et al. I l l2008h across the N90 H ii region ionized 
by NGC 602A. Thus, we do not see any evidence that the HMXB 
acquired a high space velocity following the supernova explosion. 



4 DISCUSSION 

4.1 The properties of SXP 1062 

SXP 1062 is only the third SMC X-ra y pulsar with a spin period 
larger than 1000 s jLavcock et al.l2O10l) . These slowly rotating pul- 
sars are particularly interesting because they represent a challenge 
for the theory of spin evolution of a neutron star in a close binary 
system (e.g.[lkhsanov 200^. 

According to the C orbet diagram for S MC Be-XRBs relating 
spin and orbital period l lCorbet et al.ll200f^ . we expect SXP 1062 
to have a binary period of ~300 days. Such a l ong orbital perio d 
is also expected from the /'o,.(,-EW(Hq') diagram l lReigetalJl"997h . 
from which our measured EWQla) suggests an orbital period of 
~100 days. Note however that the maximum EW(Hq') of the sys- 
tem, which probes the maximum size of the disc (and indirectly the 
orbital period), could be higher than our instantaneous measure- 
ment, so this period estimate should be taken as a lower limit. 



SXP 1062 shar es many ch aracteristics with the class of persis- 
tent Be-XRBs (e.g. lReig|l20TTh : a relatively low X-ray luminosity 
~ 10^"* - 10^^' erg s"', a slowly rotating pulsar with P > 200 s, a 
relatively flat light curve with sporadic and unpredicted increases 
in intensity by less than an order of magnitude, and a lack of iron 
lines at ~6.4keV indicative of small amounts of material in the 
vicinity of the neutron star. A thermal excess of blackbody type, 
with a high temperature {kT > 1 keV) and a small emission area 
(R < 0.5 km, consistent with a hot spot at the polar cap of the neu- 
tron star) has rece ntly been suggested as a c ommon feature of per- 
sistent Be-XRBs l lLa Palombara et aill2009h . Such a component is 
also identified in the X-ray spectrum of SXP 1062, for which we 
infer a size of ~ 0.4 km for the blackbody source from our best-fit 
models (section |2Tt . 

4.2 The environment of SXP 1062 

A shell nebula is detected around S XP 1062 in the M agellanic 
Cloud Emission-line Survey (MCELS; [Smith et al.|[T99^ Ha im- 
age a nd in the higher resolution CTIO 4m MOSAIC jMuller et"ail 
11998") Ho- image (see Figure |5}. This shell is also detected in the 
MCELS [O m] image, but hardly in [S n]. We adopt below a dis- 
tance modulus of 18.7. The shell is not uniform, with radii varying 
from 75" at the sharp rim on the northeast to 94" at the difl'use edge 
on the southwest. The peak surface brightness of the northeast rim 
is ~ 6 X 10"" ergs cm"^ s"' arcsec"-, corresponding to an emis- 
sion measure of 30 cm"^ pc. The width of the rim suggests that the 
shell thickness (AR) is 5-10% of the shell radius (R). The longest 
emitting length at the shell rim is 2R[l - (A/?/i?)^]''^; thus, the rms 
density in the shell is 1.3±0.3 H-atom cm""*. The gas mass in the 
shell is 250±100 Mg. This large mass indicates that the shell gas 
must be dominated by interstellar material. 

The shell morphology resembles supernova remnants (SNRs) 
in the Magellanic Clouds. The X-ray images also suggest that dif- 
fuse X-ray emission possibly associated with a SNR may be present 
in the vicinity of SXP 1062 (Oskinova et al., in preparation). As the 
shell mass is much larger than the typical supernova ejecta mass, we 
assume that the shell is a SNR in the Sedov phase. The kinetic en- 
ergy in the shell would be ~30% of the supernova explosion energy. 
Adopting a canonical explosion energy of 10^' ergs, the current 
shell kinetic energy is 3 x 10^" ergs, and implies a shell expansion 
velocity of 350± 100 km s"' and an age of 0.4(i?/ V')=(2-4)x 10'' yr. 
This age is much larger than the cooling timescale of this low den- 
sity gas, justifying an adiabatic shock for the Sedov phase. The 
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Figure 5. MOSAIC Ho- image and MCELS Ho-, [S n], and [O in] images 
of a region centered on tlie position of 2dFS 3831=SXP 1062 and showing 
the shell nebula detected around the target. 



pre-shock interstellar gas density, 1/4 the shell density, is 0.3 ± 0. 1 
H-atom cm"^, consistent with the low density expected in the SMC 
Wing. The bright [O in] emission can be easily produced by a 350 
km s"' shock (cf.;Hartigan et al. 1987). The ionizing flux of the BO- 
0.5III star can easily photoionize the shell gas and its surrounding 
pre-shock medium; furthermore, the diffuse [O iii] emission to the 
northeast exterior of the shell indicates the existence of a harsh radi- 
ation field; therefore, the ionization stage of the pre-shock medium 
may be too high to produce strong post-shock [S n] emission. We 
thus conclude that the shell nebula detected around SXP 1062 is 
most likely a SNR. The extent of the shell, which has not reached 
the NGC 602 region, and its young kinematic age suggest that this 
supernova event did not trigger the formation of NGC 602. 



5 CONCLUSION 

We reported the discovery of a Be-XRB, SXP 1062, containing one 
of the slowest rotating X-ray pulsars in the SMC. We suggest that 
the shell nebula detected around this object is a SNR. Its estimated 
kinematic age of (2 - 4) x lO'' yr is probably the age of the pulsar. 
This is, to our knowledge, the first discovery of a pulsar associated 
with a supernova remnant in the SMC. 
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